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Abstract 

We give a sufficient condition on a family of radial parametrized long-range potentials for a 
compact local minimality of a given d-dimensional Bravais lattice for its total energy of inter¬ 
action created by each potential. This work is widely inspired by the paper of F. Theil about 
two dimensional crystallization. 
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1 Introduction 

As explained in [2], the crystallization problem, that is to say to understand why the particles 
structures are periodic at low temperature, is difficult and still open in most cases. Theil exhibited 
in [7] a radial parametrized long-range potential with the same form as the Lennard-Jones potential 
such that the triangular lattice is the ground state of the total energy in the sense of thermody¬ 
namic limit. This kind of potential, parametrized by a real number a > 0, is bigger than a~^ close 
to the origin, corresponding to Pauli’s principle, has a well centred in 1 and a 2a width, its second 
derivative at 1 is strictly positive and its decay at infinity is r i->- ar~'^. Thus, as small is a, as 
close to 1 is the mutual distance between nearest neighbours of the ground state configuration and 
as the interactions between distant points are negligible. 

In this paper, our idea is to present a family of parametrized potential very close to this one, 
with the most natural possible assumptions, such that a given Bravais lattice L of M'’* is a N- 
compact local minimum for the total energy of interaction. This kind of local minimality is called 
“A-compact” because, given a maximal number N of points that we want to move a little bit, 
there exists a maximal perturbation of the points which gives a larger total energy of interaction, 
in the sense that the difference of energies is positive. Moreover, as small is the parameter, as 
large the number N can be chosen. We are strongly inspired by Theil’s potential, keeping only 
local assumptions and strong parametrized decay. Eurthermore, our work can be related to that of 
Torquato et al. about targeted self-assembly mu where they search for radial potentials such that 
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a given configuration - more precisely a part of a lattice - is a ground state for the total energy of 
interaction. 

The aim of this paper is to give a generic construction of the family of potentials without tak¬ 
ing into consideration the specific symmetries of L or a particular pertubation of points. 

After defining the concepts and our parametrized potentials, we give the theorem, its proof and 
some important remarks and applications. 

2 Preliminaries : Bravais lattice and N-compact local minimality 

Definition 2.1. Let d G N*, (ui,...,ttrf) he a basis of and L = Zttj C 6e a Bravais 
lattice. For any X > 0, we define m(X) := tj{L n {||x|| = A}} where ||.|| denote the Euclidean norm 
and (jA is the cardinal of set A. Moreover, we call Ai := min{||x||;x G L*}, where L* = T\{0} and 
X 2 ■= min{||x||; ||x|| > Ai,x G L}. 

Furthermore, for a Bravais lattice L and n > d, we define the following lattice sums: 

Clin) — Ikir*" and Chin) — ^ (||x|| - Ai)"”. 

x^L x^L 

Definition 2.2. Let L C be a Bravais lattice, B C L a finite subset and a he a real number 

such that 0 < a < Ai/2. We say that i?" is an a-compact perturbation of B if 

\tb G B, 3!6“ G such that \\b - 6“|| < a. 

Moreover, if B°‘ is an a-compact perturbation of B C L, we write L^{B) := {L\B) U B°‘ the 
perturbed lattice. 

Definition 2.3. Let d G N*. We say that V : M is a d-admissible potential ifV is a 

function and, for any Bravais lattice L C 

E l''(ll^■ll)l + E Il‘'llir'(l|i1l)l + E ll‘'ltlr"(l|a>ll)l < +~- 

x<=L* x<=L* x<=L* 

Remark 2.1. If, for any k G {0,1,2}, |l/(^)(r)| = 0{r~P>‘), Pk > d k, then V is d-admissible. 

Definition 2.4. Let L be a Bravais lattice o/M*^ and V a d-admissible potential. Let N G N*, we 
say that L is a N-compact local minimum for the total V-energy if for any subset B C L 
such that ^B < N, there exists oq > 0 such that for any a G [0, oq) o,nd any a-compact perturbation 
ofB, 

AUV;B):= ^ ^ y(||6« _ y||) _ ^ ^ y(||6 _ 3;||) > 0. 

yeL°^(B) b£B 

x^b 

3 Pressure, parametrized potentials and main result 

Definition 3.1. Let V be a d-admissible potential and L CW^ be a Bravais lattice, then we define 

the pressure of L submitted to V by : 

V{L,V) :=-Yl ll^ll^'dl^lD- 

x&L* 
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Remark 3.1. Actually, as explained for instance in [T], if Ey[L\ := i/(||x||) is the energy per 

x^L* 

particle of L submitted to V, i.e. the free energy at zero temperature, then, by usual thermody¬ 
namics formula, we define pressure P{L,V) by 

P(L,V) = WV'(W), 

xGL* 

where A is the area per particle of L, that is to say V{L, V) = 2AP{L, V). 

Definition 3.2. Let L C be a Bravais lattice. We call (L, 6)—family every set of d-admissible 
(L, 6)—potentials Vg : ^ M, indexed by 6 € [0, 6o], 0 < 9o < Ai/2, and satisfying the following 

conditions : 

1. Small pressure condition: there exist positive real numbers C,y such that for any 9, 

\V{L,Ve)\<C9^+^^; 

2. Parametrized fast decay : 3ro G [Ai, A 2 ), 3e > 0, 3p > d such that for any r > ro, 

\Vf{r)\ < 

3. Local convexity around first neighbours : it holds Vf{r) > rj > 0 on a neighborhood of 
the first distance Ai, uniformly on 9. 


THEOREM 3.2. Let L C he a Bravais lattice, then for any N gN*, there exists 0o > 0 such 
that for every {L, 9)—family {Vg)g<gQ and every 9 G [0, 6 * 0 ], L is a N-compact local minimum for the 
total Vg- energy. Furthermore, in this case, the maximal perturbation ao can be chosen equal to 9. 

Proof. Let L be a Bravais lattice of Let iV G N* and B := { 61 ,..., 67 V} <G L. Let oq be such that 
0 < oq < Ai/2 and be an OQ-compact perturbation of B. For any 1 < i < A^, 

for any y G L°‘°{B), y 7 ^ and x G L such that \\x — y\\ < oq, we define 

ai,x ■■= lie -y|| - \\bi-x\\. 

Obviously, we have, for any 1 < i < A", for any 0 < oq < Ai/2 and any x G L, 

\ai,x\ < 2q;o- (3.1) 

We assume, without loss of generality, that max \ai^x\ = 2 a 0 ) left to decrease ao- We set 9 G [0, Ai/2) 

i,x ’ 

and Vg a (L, 0)—potential. We have 

N N 

Al°{Vg;B) = Y, E ^^^(lie- 2 /ll)-EE^^(ll^*-^ll)- 

1=1 yeL'^O(B) i=l 

By Taylor expansion, we get, for any 1 < i < N, for any y G L°‘°{B), y fix 6 “°^ and x G L such that 

Ik - 2/11 < ao, 


VeiWbr-: 


>Vg{ 


— X 


+ Cli,xVg{'fi)i — x||) -|- 


aj 


-W" 
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for suitable S {\\bi — x|| — \ai^x\, ||^j — x\\ + |ai,a;|)- Hence we obtain 

N N 

Al°iVe-, - a^ll) + 2 E E 

i=l XGL i=l x£L 

x^bj^ x^bj^ 

We split interactions into two parts : the short range and the long range. For any 1 < i < A^, we 
set 

S'l := {x G L] ||x — bi\\ = Ai} and C\ := {x € L; ||x — 6j|| > Ai}. 

Furthermore, as we assume that for all r > tq \V 0 {r)\ < and Vq goes to 0 at +oo, 

because it is a d-admissible potential, we have, by a simple argument, that |V^^(r)| < 02 +£j,-p-i £qj, 
all r > tq. 

As L is a Bravais lattice, it holds, for any z, '^eiWbi ~ — ic|| = —V{L,Vg). Therefore, 

x£L\{bi} 


N 

Y W,xF 0 (|| 6 i -x||) 

i=l 


N 


u'(Ai) E E II''?" - i'll 


, *=i xesi 


N 


+ E E VM,-x\\)\\br-y\\+NViL,Ve). 

*=1 x£Ci 


N 

We remark that, writing Sq : = EE \\bf° — y||, and by dehnition of V{L, Vg), 

*=i xeSi 


F,'(Ai)So + NV{L, Vg) = ViL, Vg) {N - Eo(Aim(Ai))-^) 


( 


x^L 

|a:||>Al 


Eo(Aim(Ai)) ^ 


By (j3.ip . we have, for any 1 < i < A”, for any x € Ai — 2ao < \\bf° — y\\ < Ai + 2 q!o, and we get 

m(Ai)A(Ai — 2 q;o) < Eq < m(Ai)A(Ai + 2ao)- 

Thus, we show that V{L,Vg) (N — Eo(Aim(Ai))“^) > — 1 -El-2. > —by as- 

Ai Ai 

sumption. Therefore, we obtain 

H,'(Ai)Eo + NV{L, Vg) > -^A0i+^ao - . 

Hence, for first order terms, we get 

E E v^m - x||) > -^Nd^+>^ao - 2CUp)Ne^+^ - 2 + (Up + 1)] Ne^+^ao. 

i=l xGL 1 \ 1 / 

x^bi 


N 

For the second order terms, as max \aix\ = 2 q; 0 ) we have afx > 4aQ, and we obtain 

i,x ’ G—^ ^^ ’ 

*=1 x&S]^ 


1 

2 


N 

E E “A yg'Yx) > 2r/ag - 29'^^^alNC, l{p + 2). 

i=\ xgl 

x^bi 
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Finally we get, for any 0 < ao < Ai/2, 

Al°{Ve-,B) >27jal - 26^+^alNCL{p + 2) - —NO^+^^ao - 2a(p)iV02+^ 

M 

-2(^^ + C2(p+1)) NO^+^ao 
= 2r]al - N {Ae^+^ + BO^+f^ao + CO^+^ao + D9^+^al) , 
where positive real numbers A, B, C, D depend only on L. Given 6 G [0, Ai/2), if oq = 9, then 

A^iVg; B) > 27/02 _ ^ (^02+e ^ ^^2+/. ^ (jg-i+e ^ ^^4+6^ _ ^3 2) 

As 7/ > 0, there exists 9q € [0, Ai/2), depending on N, sufficiently small such that for any 9 € [0, 0o] 
and for any a € [0,0], A.J^{Vq]B) > 0 and then L is a A^-compact local minimum for the total 
V^-energy, for any (L, 0)—potential of the (L,0)—family {ye)e<eo- D 

4 Remarks 

1. Isothermal compressibility. It is usual (see Bi) to define the isothermal compressibility 
kt, from pressure P (see Remark 13.ip . by 

= n E [ii-iPv'^'dWD - ii-id^^dwi)]. 

x&L* 

where A is the area per particle of L and [L] its energy per point. We know that kt > 0 (see [SI 
Section 5.1]). Actually, that follows here from assumptions on Vg, if 0 is sufficiently small. Indeed, 
by assumptions, we have, for 0 sufficiently small, 

— = E [ll»:lt'7;'d|i|l) - lkl|Vi'{||i||)] > \\m{\gn - «"+'G(!>) - CO'+I* > 0. 

it- 

2. Zero pressure condition and local minimality among dilated of L. Let us assume here 
that, for any 0 G [0,Ai/2), 

R(L,0):=-^||x|lR,'(11x11) = 0, (4.1) 

xgL* 

which is thermodynamically natural at zero temperature, for instance if L is the cooling of an ideal 
gas. Now if we consider f : r ^ Eyg[rL\, we get, by ()4.1I) and kt > 0 (see previous remark), 
/'(I) = 0 and /"(I) > 0, i.e. L is a local minimum of L i—>• Evg[L] among its dilated, which 
seems natural if L is a A^-compact local minimum for the total V^-energy for N arbitrarily large, 
and this is actually assumed in Theil’s paper [Tj. 

However, the reverse is false. A Bravais lattice can be a local minimum among its dilated for 
the energy per point but not a A^-compact local minimum for the total energy. For instance, if 
d = 1, L = Z, Ai = 1 and 1/ defined by V{r) = 0 for r > 5/2, V'{1) = V'{2) = 0, V"{1) = -1 
and V''{2) = 1/3. We have ^ |x]lL'(|x]) = 0 and ^ |x]2r"(|x]) = 2/3 > 0, hence Z is a local 

minimum among lattices of the R-energy per point. For a > 0, we estimate, by Taylor expansion, 

A“(R;L)= ^ [R(|x-a])-R(]x|)] 

= a^V''{l) + a‘^V"{2) + a‘^(P{a) = a^{-2/2, + (/.(a)). 
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where 0(a) goes to 0 as a —>■ 0. Hence for a < ao sufficiently small, —2/3 + 0(a) < 0 and Z is not 
a 1-compact local minimum of the total H-energy. 

3. Effects of parameters e,fi,p and t]. By (13.2p . our assumptions on V 0 give indications 
about the stability of lattice L: 

• Range : a large p, C, e 01 p allow to take a large perturbation ao for fixed N, i.e. a better 
decay at infinity implies a stronger stability of the lattice; 

• Second derivative around nearest-neighbours distance : a large rj also allows a large pertur¬ 
bation ao for fixed N. Typically, a narrow well around Ai “catches” the first neighbours of 
the minimizing configuration at distance Ai. 

4. Difference between the decay after the first distance and the pertnrbation. We can 

see that 6 ^'^^ « 60 , i.e. the maximum decay is really smaller than the maximum perturbation and 
this allows not to assume a local behaviour of Vg around Ai with respect to 9, as in Theil’s work. 
Obviously, if 0 = 0 then Vb(r) = 0 for any r > vq and V)((Ai) = 0, therefore Ai is a local minimum 
of Vq and the potential is short-range : only the first neighbours interact and the A^-compact local 
minimality is clear for any N with a perturbation ao as small as N is large. 

5. A kind of Cauchy-Born rule. Our result can be viewed like a justification of a kind of 
Cauchy-Born rule (see mm)- Indeed, if we consider a solid as a Bravais lattice L where the inside 
is a finite part of L with cardinal N and the rest is its boundary, a small perturbation of the inside, 
depending on N, increases the total energy of interaction in the solid. That is to say that the inside 
of the solid follows its boundary to a stable configuration. 
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